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ABSTRACT 


An  efficient  numerical  procedure  for  the  transient  dynamic  analysis 
of  elasto-plastic  shells  is  introduced.  A simple  shell  is  analyzed,  and 
the  results  achieved  are  compared  against  an  existing  code. 
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NOMENCLATURE 


Area  of  element  I 

[B]  Strain-displacement  matrix 

c^  Distance  between  neutral  axis  of  a stiffener  and 

the  reference  surface  of  the  shell. 


[D] 


Tangent  moduli  matrix 


(e> 


Strain  vector 


{e^},{e2),{e^},{e^} 

Cxx’  eyy’  6xy 
ell*  e22’  e12 


Strain  vector  in  element  1 

Strain  vectors  in  the  4 regions  within  an  element; 
components  of  vector  {e^}. 

Components  of  strain  in  cartesian  coordinates; 
components  of  vector  {e^},  {e2},  etc. 

Components  of  strain  in  an  orthogonal  coordinate 
system. 


F’  V fl 


F , F 
s’  M 


IM’  XNM 

k , k 
yy  xy 


kll’  k22’  k12 


Subsequent,  initial,  and  limit  yield  functions 

Absolute  values  of  the  gradients  of  the  yield 
function  F. 

Force  vector  of  element  i 
Metric  coefficients 
Stress  resultant  invariants 

Components  of  curvature  within  an  element,  in 
cartesian  coordinates;  components  of  the  strain 
vector  (e^},  (e2)  etc. 

Components  of  curvature  in  an  orthogonal  coordinate 
system. 


Mll’  M22’  M12 


M 


ij 


Distance  between  stiffeners 

Moment  per  unit  length;  components  of  the  stress- 
resultant  vector 

Strain  hardening  parameters 

Total  number  of  elements  with  the  structure. 


Nll’  N22’  N12 


Normal  force  per  unit  length;  components  of  the 
stress-resultant  vector. 


✓ 


f 


J 


{P}  Vector  of  external  forces  acting  on  the  nodes 

{ p } Surface  loading  per  unit  area 


{s> 

M 


w 


wr  w2  ’ ’ *wi2 


x»  y 


Nodal  displacement  vector  of  element  i. 

Radii  of  curvature  along  the  principal  directions 
of  the  shell  in  cartesian  coordinates. 

Radii  of  curvature  along  the  principal  directions 
of  the  shell  in  orthogonal  coordinates. 

Stress-resultant  vector 

Displacements  of  a point  within  an  element. 

Tangential  displacements;  components  of  the 
vector  {u}  in  section  II. 

Tangential  displacements  in  the  x direction  of 
the  4 nodes  contiguous  to  an  element;  components 
of  vector  {q}  in  section  IV. 

Velocity  vector  of  node  j at  time  step  n. 

Tangential  displacements  in  the  y direction  of  4 
nodes  contiguous  to  an  element;  component  of  vector 
{q}  in  section  IV. 

Normal  displacements;  components  of  vector  {u}  in 
section  II. 

Normal  displacements  of  the  12  nodes  around  an  element; 
components  of  vector  {q}  in  section  IV. 

Cartesian  coordinates  x = hj^,  y = ^2^2 


At 

Ax,  Ay 
5e 

6u 


^1*  ^2 
P 

4>.  ‘f’jt  <t>2 


Time  step 

Distance  between  nodes  in  cartesian  coordinates. 
Increment  in  strain 
Increment  in  displacement 

Orthogonal  coordinates  along  the  principal  curvatures. 
Mass  density  per  unit  area  of  shell  surface. 

Twist,  rotations  with  respect  to  and  • 
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I INTRODUCTION 

This  report  is  part  of  a combined  theoretical-experimental  study, 
under  a joint  DNA/NAVSEA/ONR  program,  which  is  aimed  at  evaluating  and 
increasing  submarine  hardness  to  underwater  explosions.  Previous  reports 
in  the  series  were  aimed  primarily  at  studying  shock  effects,  in  the 
elastic  range,  on  internal  equipment.  The  present  report  is  part  of  an 
effort  to  analyze  submarine  lethality  problems,  where  large  elasto- 
plastic  motions  of  the  hull  may  occur  under  the  action  of  long  duration 
full  envelopment  shock  loadings,  which  are  produced  by  nuclear  explosions. 

The  first  report  in  the  lethality  series.  Ref.  [1],  involves  the 
formulation  of  an  elasto-plastic  theory  for  stiffened  shells.  The  present 
report  describes  a numerical  procedure  for  the  dynamic  analysis  of  elasto- 
plastic  shells  using  a direct  integration  in  time.  While  the  examples 
presented  here  are  for  shells  in  vacuo,  the  Doubly  Asymptotic  Approximations 
(DAA)  is  currently  being  implemented  into  the  code  in  order  to  treat  the 
fluid  interaction  problem. 

The  objectives  of  this  work  are  the  following: 

1)  A realistic  modeling  of  complex  shells  of  arbitrary  geometry, 
including  stiffeners  and  internal  structure. 

2)  The  analysis  of  shock  phenomena  with  high  frequency  components 
in  their  spectrum. 

3)  The  modeling  of  elasto-plastic  material  behavior. 

4)  The  capability  of  taking  into  account  large  displacement  gradients 
in  order  to  analyze  dynamic  buckling  conditions  and  post-buckling 
behavior. 

The  realistic  modeling  of  complex  structures  for  shock  phenomena 


results  in  very  large  problems,  with  many  degrees  of  freedom.  This  fact 
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and  the  nonlinearities  implied  by  the  conditions  (3)  and  (4)  point  to 
a finite  element  or  finite  difference  method,  together  with  the  direct 
integration  in  time  of  the  equations  of  motion,  as  the  most  promising 
approach.  The  nonlinearities  of  the  problem  makes  the  normal-mode 
expansion  method  unsuitable  for  the  present  purposes. 

A survey  of  the  existing  numerical  methods  of  analyzing  elasto- 
plastic  shells  for  shock  problems,  undertaken  at  the  beginning  of  this 
work,  revealed  that  additional  development  was  needed.  There  are 
several  aspects  to  the  existing  methods  which  make  them  inappropriate 
in  meeting  the  analysis  requirements  of  the  present  work.  Specifically: 

(a)  A refined  element,  employing  high-order  approximating  functions,  has 
proven  to  be  an  efficient  approach  to  the  static  analysis  of  shells 
and  the  determination  of  low  natural  modes  and  frequencies.  The 
accuracy  of  such  an  element  allows  for  the  use  of  larger  elements, 
reducing  the  total  number  of  degrees  of  freedom  required  for  the 
entire  shell.  In  contrast  to  this,  in  shock  and  wave  propagation 
problems,  a large  number  of  mass  points  is  the  only  way  of  achieving 
an  accurate  solution,  even  if  the  element  itself  is  relatively  simple 
It  appears  that  the  shell  elements  with  "condensed"  stiffness 
matrices,  i.e.  with  massless  and  loadless  internal  nodes  eliminated 
by  a static  condensation  procedure,  are  especially  unsuitable  for 
this  purpose. 

(b)  The  masses  associated  with  the  rotational  degrees  of  freedom  present 
in  the  finite  element  method  are  usually  very  small.  This  causes 
numerical  difficulties  in  an  explicit  scheme.  If  the  "rotational 
masses"  are  neglected,  the  solution  consists  of  the  integration  in 
time  of  the  translational  degrees  of  freedom,  together  with  the 
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solution  of  a system  of  "static"  equations  for  the  rotational 
degrees  of  freedom.  Since  these  equations  are  nonlinear,  their 
solution  (which  must  be  repeated  at  each  step  of  the  integration 
in  time)  would  require  an  enormous  increase  in  the  number  of 
computations. 

(c)  For  problems  involving  a structure  subject  to  cyclic  loads  in 
the  plastic  range,  the  determination  of  the  stress  resultants 
in  existing  methods  is  accomplished  by  first  computing  the 
stress  components  at  several  locations  through  the  thickenss 
of  the  shell,  and  then  computing  the  stress  resultants  by 
numerical  integration  of  the  stresses  through  the  thickness  of 
the  shell.  While  the  computations  involved  are  simple,  the 
amount  of  stored  data  becomes  prohibitively  large  for  structures 
involving  many  elements. 

The  approach  adopted  in  the  present  work  is  believed  to  be  free 
from  the  above  mentioned  drawbacks.  The  main  aspects  of  this  approach 
and  their  current  status  are  discussed  in  the  following  sections. 
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II  BASIC  SHELL  EQUATIONS 


The  kinematic  equations  of  the  shell  theory  employed  in  this  work 
correspond  to  the  Donnell-Vlasov  nonlinear  theory,  with  the  option  open 
of  refining  them  to  achieve  Sander's  theory.  In  orthogonal  coordinates 
along  the  principal  curvatures,  the  strain-displacement  relations  read 

9ul  1 9hl  ._w,l,21,2 

G11  = h^  hxh2  U2  R!  2 *1  2 * 


du2  . 1 8h2  . w . 1 , 2 . 1 ,2 

e22  = 1^51 r2  + W;  U1  + + 2 *2  + 2 * 


3u2  3ux  l 3h1  x 9h2 

2e12  " h^  + h23C2  " hxh2  3^2  U1  " h2h2  3^  U2  + *1*2 


k - 34>1  + 1 ahz 

11  “ h^  h^  3^  <J>2 


k ^2  , ! 8h2, 

*22  h23C2  + hLh2  352  * 1 


2k  = 34>2  + 301  1 9hl  * _i_  +(I  _!_) 

12  h13C1  h23C2  ' hxh2  8C2  V1  “ hxh2  <P2  + 


* _ 3w  . _i 

*1  h13^1 


3w  _2 

h23C2  R2 


• 1 1 r 3 ,,  . 3 ,,  .. 

* ‘ 2 hxh2  ^ (h2u2)  ■ 3^2  (hlul)] 
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The  underlined  terms  represent  geometric  nonlinearities;  the  terms 
with  broken  underlines  are  those  of  Sander's  theory. 


The  stress  resultants  are  defined  as 

h/2 


Nid  = 


-h/2 

h/2 


Sij  dZ 


*ij  = J Sij  2 dZ 

-h/2 


with  i = 1,  2;  j = 1,  2. 

The  theory  is  completed  by  writing  the  equations  of  equilibrium 
which,  for  the  present  purposes,  are  assumed  in  the  form  of  the  Principle 
of  Virtual  Work.  With  the  notation 


{u}  = 

(v 

u2,  w)T 

{s}  = 

(Nn, 

N22’  N12’  MH*  M22’  M125 

{e}  = 

(eil’ 

e22’  2e12’  kll’  k22’  2kl2) 

{p}  = 

(Pi* 

^T 

P2.  P3) 

standing  for  the  surface  loading.  The  Principle  of  Virtual  Work  reads 


{s}^{6e}  dS  — J {p}^{6u}  dS  + [ p{u}^{6u}  dS  = 0 

JS  JS  JS 


where  p is  the  mass  density  per  unit  area  of  the  shell. 
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III  SHELL  CONSTITUTIVE  EQUATIONS 

The  shell  constitutive  equations  are  the  relations  between  the 
rates  of  the  stress  resultants  and  the  rates  of  the  shell  strains.  In 
the  matrix  notation  introduced  in  the  preceeding  section,  they  are 

U}  = [D]  { e } (3.1) 

where  [D]  is  the  so-called  elasto-plastic  tangent  stiffness  matrix. 

The  explicit  from  of  the  above  relation  is  based  on  an  elasto- 
plastic  theory  for  shells  presented  in  a separate  report  (Ref.  [1]). 

The  aforementioned  theory  is  analogous  to  classical  theories  of 
plasticity;  consisting  of  a yield  condition,  a strain  hardening  rule, 
and  a flow  rule.  It  differs  from  classical  elasto-plastic  theories  in 
that  the  yield  surface  is  defined  in  terms  of  the  stress  resultants  (s) 
instead  of  the  stresses.  This  avoids  the  necessity  of  computing  and 
storing  stresses  through  the  thickness  of  the  shell. 


The  stress  components  at  the  top  and  bottom  surfaces  of  a solid 
shell  can  be  expressed  in  terms  of  the  stress  resultants  as: 


a 


ij 


where  h is  the  thickness  of  the  shell,  and  the  plus  sign  applies  to 
the  top  and  the  minus  sign  to  the  bottom  of  the  shell. 

An  expression  for  the  initial  yield  surface  is  constructed  by 


(3.2) 


resulting  in 


F0  = XN  + XM  ± 2 rNM  " 1 


where 


and 


TN  5 ~2<N11  + N22  ' NU  N22  + 3N12> 

No 

V E 72<"ll  + «22  - H11  «22  + 3H12> 

Mo 

*NM  = NqM0(N11  MH  + N22M22  " 2 N11N22  “ 2 N22M11  + 3N12M12) 


N0  = °Oh  * M0  = °Oh  /6 


A limit  surface  is  assumed  which  is  also  a linear  combination  of 
Ijj*  Ijj,  and  I . Coefficients  for  the  three  terms  are  determined 
empirically  in  order  to  produce  a good  approximation  for  the  limit  surface. 
The  expression, 

FL  = IN  + 9 ^ - 3 yTT  INM 

which  represents  the  limit  condition  exactly  for  the  three  extreme  cases  of 
(1)  only  membrane  forces,  (2)  only  bending  moments  and  (3)  = N22  with 

= M22>  was  chosen. 

A hardening  rule  is  generated  in  the  following  manner.  A variable 
yield  surface  of  the  form 


F = I +I*+ai  =1 
N M - NM  i 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 
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is  assumed  where 


Hi  “ 2 [M11  " Mll)  + (M22  " M22J 
M0 

(Mll  _ Mll^M22  ” M22^  + 3^M12  “ M12^  ^ 


The  quantities  M^,  which  represent  "hardening  parameters",  are  defined 
by  the  following: 

If  F-i  and  ^ ^ =•  0: 

Mo  f82  " 

“u  - 2(1  - v 4 f2  % 

M 

I£F<1  or  ^ *lj  + 3^  "ij  i 0: 
d M^*  « 0 


The  symbols  F and  F__  are  defined  as 

S M 


'o  aN^ 

+ <"0  3H22» 

2 

3F  . 

> ) 

+ <M„  J?-) 

+ (M„ 

2 1/2 


2 2 1/2 

FM=  [(M0  S£>  + (M0^>  + % 9S^>  1 


Fg  is  the  absolute  value  of  the  vector  grad  F,  in  a dimensionless 

formulation;  F is  the  part  of  grad  F which  corresponds  to  the  bending 
M 

moments  only. 


A more  detailed  discussion  of  the  theory,  along  with  some  numerical 
results,  are  presented  in  Ref.  [1]. 


I 
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IV  DISCRETIZATION 

The  discretization  procedure  adopted  has  two  aspects  in  common  with 
the  classical  finite  element  method.  First,  a variational  derivation  is 
used  for  the  equations  of  motion.  Second,  a quadrilateral  element  is 
used,  in  which  the  stresses  and  strains  are  computed  from  the  nodal  dis- 
placements of  nodes  surrounding  the  element. 

A quadrilateral  shell  element  defined  by  four  corner  nodes,  with 
each  node  having  three  translational  and  no  rotational  terms,  does  not  have 
enough  degrees  of  freedom  to  represent  bending  behavior  (second  derivative 
terms).  In  order  to  have  the  required  additional  degrees  of  freedom, 
eight  nodes  not  contiguous  with  the  element  are  also  used  in  computing 
bending  terms  (Fig.  1). 


Each  element  accesses  twelve  nodes  and  has  twenty  degrees  of  freedom; 
three  translational  degrees  of  freedom  for  each  of  the  four  inner  nodes, 
and  one  degree  of  freedom  (displacement  normal  to  the  surface)  for  each  of 
the  eight  exterior  nodes.  Thus,  the  nodal  displacement  vector  of  an 
element  i is 


(V  u2,  u3,  u4,  vr 


(4.2) 


For  the  purpose  of  computing  the  stiffness  of  an  element,  the  area 
of  the  element  is  divided  into  four  regions  (Fig.  2).  One  set  of  strains 
is  computed  for  each  of  the  four  regions.  The  strain  in  element  i is 


{e}^  = ({e^,  (e2),  {e^},  {e^})  wherefe^}  is  the  strain  in  region  1,  and 

T 


(e  ,e  , 2e  , k.  , k , 2k  ) 
xx  yy  xy  xx  yy  xy 


At  the  present  time,  in  order  to  analyze  a series  of  simple  check 
problems,  the  discretization  has  been  limited  to  the  case  of  small  strains 
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expressed  in  cartesian  coordinates.  In  this  simplest  version  of  shell 
theory,  the  strains  have  the  following  form: 


3u  . w 3v  . w lr3u  . 3v. 

exx  = 3x  Rx  ’ eyy  “ 3y  Ry  ’ ®xy  “ 2f3y  3x] 


3^w  . = 3j#  . = 32w 

3x2  ’ yy  " 3y2  ’ xy  ^ 


where  R and  R are  the  principal  radii  of  curvature  of  the  shell  surface, 
x y 

The  discrete  approximations  to  c above  expressions,  in  region  1 
of  A^,  are  (for  equally  spaced  nodes) 

3u  = V-L  3v  = V!l  3u  = Vu_l  3v  = V^l 
3x  Ax  * 3y  Ay  ’ 3y  Ay  ’ 3x  Ax 


(4.2) 


J2  w„-2w  -w_  »2  w.  -2w.+w. 

3w  215  3w  416 


(4.3) 


.2  w^-w.  - w.+w, 

3 w _ _3 2 4 1 

3x3y  AxAy 


Analogous  expressions  are  formed  for  regions  2,  3 and  4.  This 
process  defines  a strain-displacement  matrix  [B]^  such  that 


{e>i  = [B] ± {q>1 


(4.4) 


The  expression  for  derivatives  is  basically  a staggered  finite 
difference  scheme,  with  first  derivatives  computed  between  nodes,  and 
second  derivatives  at  nodes.  One  should  note  that  there  is  no  displacement 


function  assumed  over  the  element  area  for  the  determination  of  derivatives, 


as  one  finds  in  a classical  finite  element  formulation. 
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This  type  of  element,  which  uses  non-contlguous  nodes,  results  In 
an  overlapping  of  adjacent  elements.  This  overlapping  produces  compli- 


catlons  at  boundaries  similar  to  those  encountered  In  finite  difference 


formulations,  but  are  absent  from  finite  element  formulations.  To 
minimize  the  difficulties  associated  with*  boundaries,  the  structure  is 
conceptually  divided  into  sheets.  Each  sheet  is  a curved  section  of 
shell  with  an  arbitrary  number  of  nodes  and  elements  (Fig.  3 ).  The 
shape  of  the  sheet  is  limited  to  a surface  that  can  be  described  by  a 
smooth  continuous  function  without  any  interior  discontinuities  in  its  slope. 
The  elements  within  a sheet  are  limited  to  a purely  rectangular  organi- 
zation (i.e.  exactly  four  elements  must  be  connected  to  each  interior  node). 


Thus  a cylinder  with  end  caps  would  consist  of  three  sheets;  a 
circular  cylindrical  sheet,  and  a sheet  for  each  end  cap  (Fig.  4) . Three 
sheets  are  required  to  specify  this  structure  because  of  the  edge  that 
occurs  between  the  cylinder  and  the  end  caps. 

By  organizing  the  structure  into  sheets,  all  the  difficulties 
associated  with  boundaries  are  isolated  into  a set  of  artificial  nodes 
around  the  perimeter  of  the  sheet.  When  several  sheets  are  connected 
at  an  edge,  the  following  conditions  are  imposed  on  the  nodes  along  a 
shared  edge: 

1)  Comparability  of  displacements  of  nodes  along  the  edge. 

2)  Comparability  of  rotations  of  the  nodes  along  the  edge. 


3)  Equilibrium  of  moments  at  the  nodes  along  an  edge. 


V EQUATIONS  OF  MOTION 


With  the  shell  surface  divided  into  N elements  A^,  the  Principle  of 
Virtual  Work  from  Section  2 becomes 


i=N 

l l 

i-1 


{s}i{6e>i  dA  - 


{p}  {6u}  dA  + 


p{ii}  (6u)  dA  = 0 


(5.1) 


If 

of  nodal 


the  variations  of  strain  are  expressed  in  terms  of  the  variation 
displacements 


{6e>1  = [B]i(6q}1 


(5.2) 


and  with 


{s}^  [B]1  dA  « {F)i 


l f {p }T{ 6u } dA  = {p}T{6q} 

i=l  [ 


i=N 

l 

i-1 


p(ii}T  (6u)  dA  = ([m]  {q})T(6q} 


(5.3) 

(5.4) 

(5.5) 


the  Principle  of  Virtual  Work  results  in  the  following  system  of  ordinary 
differential  equations 

i=N 

[m]  {q}  = ~l  {F}  + {P}  (5.6) 

i=l 

In  the  above  equations,  (q)  is  the  nodal  displacement  vector  for 
the  Structure, [m]  is  the  lumped  mass  matrix,  and  {P}  represents  the 
vector  of  external  forces  acting  on  the  nodes  of  the  structure. 
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VI  INTEGRATION  IN  TIME 

An  explicit  scheme  of  integration  in  time  has  been  selected  for 
the  following  reasons: 

(a)  An  implicit  scheme  would  require  assembly  and  inversion  of 
large  matrices  whose  elements,  in  general,  change  at  each  step 
of  one  integration.  This  aspect  would  seriously  impair  the 
ability  to  deal  with  large  problems. 

(b)  The  more  favorable  stability  of  an  implicit  scheme  would  not 
be  utilized,  since  in  most  shock  loading  problems  the  size 
of  the  time  step  is  limited  by  accuracy  requirements  rather 
than  stability  conditions. 

The  intended  incorporation  of  the  method  of  shell  analysis  resulting 
from  this  work  into  a program  for  three-dimensional  dynamics  of  the 
interacting  medium  (e.g. , TRANAL)  provides  an  additional  strong  argument 
for  the  direct  and  explicit  integration  of  the  equations  of  motion. 

The  system  of  equations  for  the  nodal  displacements  of  the  structure, 
derived  in  the  preceding  section,  is  integrated  in  time  with  the  following 
central  difference  scheme.  The  velocity  of  node  j at  time  step  n+1  is 
computed  by 

{q>?+1  = fq}?  ({PL  - l {F}.) 

J J i«j  i 1 

where  element  forces  {F } _^  are  summed  over  all  elements  i framing  into  node 
j,  nij  is  the  mass  of  the  node,  and  {P}^  are  the  externally  applied  forces 
on  the  nodes. 


I 
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VII  SMALL  STIFFENERS 


The  term  "small  stiffeners"  is  applied  to  relatively  small  stiffeners 
whose  spacing  is  small  enough  to  allow  for  their  representation  as  a 
continuous  additional  layer  of  the  shell  structure.  This  assumption 
is  identical  to  the  assumption  which  leads  to  the  treatment  of  stiffened 


shells  as 


Lc  shells  in  the  theory  of  elastic  shells. 


An  element  of  the  stiffened  shell  is  conceptually  divided  into  the 
shell  sheet  and  the  stiffener  (Fig.  5).  The  middle  surface  of  the  shell  sheet 
remains  the  reference  surface  of  the  stiffened  shell.  Within  the  shell 
sheet,  the  increment  in  strain  and  the  increment  in  stress  resultants 
are  computed  in  the  usual  manner.  That  is, 


{de}  = [Bj  {dq} 


(7.1) 


{ds}  = [D]  {de} 


(7.2) 


For  stiffeners  in  the  direction,  the  strain  increment  in  the 

shell  sheet  is  referred  to  the  centroidal  axis  of  the  stiffener.  Thus, 

(dkll) stiff  = (dkll) sheet  ’ (dell)stiff  = (deil)sheet  _Cl(dkll) sheet 
where  "c^"  is  the  distance  between  the  centroidal  axis  of  the  stiffener 
and  the  reference  surface.  The  increments  in  the  stress  resultants  for  the 
stiffeners  are  computed  using  elasto-plastic  constitutive  equations  for 
a beam,  which  are  simply  a reduced  form  of  the  shell  constitutive  equations 


(7.3) 


The  total  increment  of  the  stress  resultants  for  the  stiffened  shell 


is  the  sum  of  the  stress  resultants  for  the  shell  sheet  and  the  stress 


resultants  of  the  stiffeners  referred  back  to  the  reference  surface. 
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Thus, 


"ll  ■ (dNu>,h,et  ♦ 4 <dNll>.tlff 


"11  ■ <dMll>.h«t  + 4 <"n  - 'ldV.tl« 


where  is  the  distance  between  small  stiffeners. 


(7.4) 

(7.5) 


1 
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VIII  LARGE  STIFFENERS 

For  large  stiffeners  the  element  mesh  must  be  arranged  in  such  a way 
that  a mesh  line  runs  along  a large  stiffener.  The  stiffener  is  treated 
as  a distinct  curved  beam  element  lying  along  the  mesh  line  of  the  shell. 

The  nodal  displacement  vector  for  a beam  element  is 

{q}j  = (u1,  u2,  vx>  v2,  w1,  w2,  w3,  w4)  (8.1) 

where  u^-w^  are  shown  in  Fig.  6 . Beam  strains  are  computed  from  nodal 
displacements  using  a reduced  set  of  relations  consistent  with  those  used 
in  the  shell  (Eq.  4.3,  4.4). 

From  this  point,  the  beam  is  treated  in  exactly  the  same  manner 
as  the  small  stiffener.  That  is,  the  strains  are  referred  to  the  centroid 
of  the  beam  (Eq.7.3).  The  stress  resultants  are  computed  from  the 
strains,  and  these  resultants  are  then  referred  back  to  the 
middle  surface  of  the  shell.  The  beam  contribution  to  the  nodal  forces 
along  {qh  follows  from  the  principle  of  virtual  work. 

K-miv*  <s-2) 

s.  3 

i 


i 


IX 


IMPLEMENTATION 


The  formulation  summarized  in  the  preceeding  sections  is  in  the 
process  of  being  incorporated  into  a comptuer  code  called  SHELPLAS.  The 
present  version  of  SHELPLAS  is  limited  to  the  small  displacement  analysis 
of  shells  modeled  by  one  sheet  of  rectangular  elements.  Stiffeners  and 
elasto-plastic  material  behavior  have  been  incorporated  in  this  version. 

In  order  to  verify  the  overall  approach  incorporated  in  this  analysis, 
some  test  problems  were  run  and  the  results  obtained  were  compared  against 
results  computed  with  DYNAPLAS  (Ref.  [2]). 

A circular  cylindrical  shell  with  a radius  of  16.8125"  was  modeled 
as  shown  in  Fig.  7.  Forty  five  elements  were  used  along  the  length  of  the 
shell.  Small  stiffeners,  spaced  at  5.625  inches,  were  modeled  as  an 
orthotropic  shell.  The  large  stiffeners,  spaced  at  50.625  inches,  were 
modeled  as  a single  element.  The  left  end  of  the  shell  was  built  in,  and 
the  right  end  was  a plane  of  symmetry. 

The  model  was  subject  to  a triangular  load  (in  time),  see  Fig.  8,  of  uniform 
pressure  (in  space).  The  peak  pressure  was  670  psi,  and  the  duration  of 
the  load  pulse  (.5  millesec.)  was  close  to  the  period  of  the  breathing 
mode  of  the  shell  (.55  millisec). 

A history  of  hoop  stresses  (Oq)  at  element  it 28  was  plotted  in  Fig.  8. 

One  can  see  that  there  was  considerable  yielding  during  the  first  cycle. 

Figure  9 is  a graph  of  the  displacements  versus  time  from  both  codes. 

Peak  displacements  agree  to  within  less  than  one  percent.  Excellent 

agreement  was  achieved  for  all  points  a sufficient  distance  from  the 

large  stiffeners,  as  can  be  seen  from  a graph  of  the  displaced  shape  (Fig.  10). 
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Figures  11  through  13  show  displacement  histories  at  the  three 
nodes  closest  to  the  large  stiffeners.  As  one  approaches  the  stiffener, 
the  gradient  of  the  moment  becomes  considerably  larger.  The  differences 
in  the  results  of  the  two  codes,  as  one  approaches  the  stiffener,  were 
due  to  the  different  approximations  of  in  the  elements  of  the  two  codes. 


The  execution  time  for  both  codes  was  less  than  ? minutes. 
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X CONCLUSIONS 

As  mentioned  in  the  introduction,  this  work  is  concerned  with  the 
inelastic  response  of  submarines  under  dynamic  loadings  generated  by  a 
nuclear  explosion.  The  complexities  involved  in  modeling  submarine 
hulls  with,  major  internal  components  require  that  careful  consideration 
be  given  toward  optimizing  computer  usage  in  the  analysis.  To  this  end, 
an  elasto-plastic  shell  theory  has  been  developed  and  reported  in  Ref.  [1], 
The  present  report  describes  a numerical  procedure  and  a direct  integration 
code  which  uses  the  shell  theory  in  the  analysis  of  stiffened  cylindrical 
shell  structures. 

The  excellent  agreement  achieved  in  comparison  with  other  codes  for 
the  analysis  of  a series  of  simple  check  problems  has  verified  the  correct- 
ness of  the  overall  approach.  For  the  eventual  analysis  of  the  very 
complex  submarine  problems,  it  is  felt  that  the  present  numerical  pro- 
cedures will  prove  much  more  efficient  than  presently  available  methods. 

The  code  presented  in  this  report  is  currently  being  augmented  to 
include  the  following  aspects: 

(1)  Fluid  - structure  interaction  using  the  Doubly  Asymptotic 
Approximation. 

(2)  Large  displacement  capability  to  study  dynamic  buckling  and 
postbuckling  behavior. 

(3)  Complex  shell  structures  with  internal  components. 

It  is  currently  planned  that  this  work  will  be  sufficiently  advanced 
to  enable  us  to  make  predictions  for  the  experiments  on  stiffened 
cylindrical  structures  which  will  be  conducted  in  the  second  half  of  1977. 
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